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ABSTRACT 

An experimental and numerical investigation is presented of a H2/N2 turbulent jet flame 
burner that has a novel vitiated coflow. The vitiated coflow emulates the recirculation region of 
most combustors, such as gas turbines or furnaces. Additionally, since the vitiated gases are 
coflowing, the burner allows for exploration of recirculation chemistry without the 
corresponding fluid mechanics of recirculation. Thus the vitiated coflow burner design 
facilitates the development of chemical kinetic combustion models without the added complexity 
of recirculation fluid mechanics. 

Scalar measurements are reported for a turbulent jet flame of H2/N2 in a coflow of 
combustion products from a lean ((j)= 0 . 25 ) H2/Air flame. The combination of laser-induced 
fluorescence, Rayleigh scattering and Raman scattering is used to obtain simultaneous 
measurements of the temperature, major species, as well as OH and NO. Laminar flame 
calculations with multi-component diffusivity are presented and do not agree well with the 
experimental results. Laminar flame calculations with equal diffusivity do agree when the 
premixing and preheating that occurs prior to flame stabilization is accounted for in the boundary 
conditions. Also presented is an exploratory pdf model that predicts the flame’s axial profiles 
fairly well, but does not accurately predict the lift-off height. 
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INTRODUCTION 


Advanced combustor designs employ the back mixing of hot combustion products with 
cool reactants in order to achieve flame stabilization and keep combustor size to a minimum. 
This important feature is absent in the widely studied simple jet flame issuing into a coflow of 
air. Attempts have been made to incorporate this recirculation with bluff body and swirl flames. 
The shortfall to these approaches is that the chemical kinetics and complex recirculation and wall 
interactions are strongly coupled. 

The coaxial jet flame in a coflow of hot combustion products (vitiated coflow) simplifies 
the fluid dynamics by emulating the preheating and premixing of recirculation without the actual 
recirculation. As a result, the decoupled chemical kinetics can be modeled with more detailed 
chemical mechanisms. The design allows for intense turbulent mixing while maintaining a 
stable flame, thus providing the opportunity to examine turbulent mixing and combustion under 
the flow conditions typical of, or beyond, advanced combustors. Most importantly, the novel, 
open configuration of the vitiated coflow burner provides both optical access and well-defined 
boundary conditions, thus making it amenable to optical diagnostics and computational 
explorations. In addition, flameless oxidation (“FLOX”) can easily be explored with this burner 
[Wiinning & Wunning 1997, Plessing et al. 1998], 

The principle objective of our research efforts is to characterize and model lifted jet 
flames of simple hydrocarbons (CFLf). As an incremental step towards that goal, the current, 
relatively simple lifted hydrogen jet flame is examined. There have been several comprehensive 
studies of hydrogen jet flames that included experimental and numerical efforts [Barlow et al. 
1994, 1996, Meier et al. 1996, Neuber et al 1998, Chen et al. 1996]. Furthermore, there are at 
least three data sets [Brockhinke et al. 2000, Tacke et al. 1998, Cheng et al. 1992] that provide 
comprehensive information for lifted H 2 flames in a coflow of air. The vitiated coflow provides 
a fairly small increment of complexity with benefits that can contribute to the understanding of 
turbulent combustion fundamentals of hydrogen and subsequently, simple hydrocarbons. 

Multiscalar measurements, even with modest accuracy, constitute a useful test for 
combustion models [Warnatz et al. 1999], This report focuses on the multiscalar measurements 
of a lifted F12/N2 jet flame in a vitiated coflow. These results are compared to laminar flame 
calculations with equal diffusivity and full transport and an exploratory PDF model. 

EXPERIMENTAL AND NUMERICAL METHODS 


The combustor consists of a central F12/N2 jet and a vitiated coflow as shown in Figure 1 . 
Note that the flames in Figure 1 are CPU flames. The coflow consists of combustion products of 
a highly turbulent, lean premixed F^/Air Flame. The coflow flame is stabilized on a perforated 
disk. In order to stabilize a highly turbulent flame, the blockage of the perforated plate should be 
80-90%. The coflow must be highly turbulent, or have enough momentum so as not to collapse. 
A detailed design package has been prepared and is accessible over the World Wide Web [Cabra 
2000]. 

The nozzle exit is high enough (7cm) above the perforated plate so that a uniform flow 
field can be assumed, based on the idea that the small-scale turbulence generated by a perforated 
plate is quickly dissipated. The temperature and oxygen level of the coflow can be adjusted; the 
operating range of the burner has previously been presented [Cabra et al. 2000], 
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Figure 1. 

Vitiated coflow combustor schematic, and photo of CFE-Air Jet in a CFU-Air Coflow. 

Experiments were conducted on one jet flame that consisted of 25% H 2 in N 2 by volume 
with a total flow rate of lOOslm. The Reynolds number of the jet was 24,000. The coflow was 
products from a H 2 /Air premixed flame ( 0=0.25 ) with a temperature of 1035K and a total flow 
rate of 2325slm. The height of the base of the lifted flame from the nozzle exit was 
approximately z/d=10 and the height of the flame was z/d=30. These characteristics are similar 
to those reported by Tacke et al. 1998 on a lifted H 2 /N 2 jet of the same composition (25%/75%) 
with a similar Reynolds number of 17,000 but with a cool, slow moving coflow (0.2 m/s). More 
information on the geometry and the flow conditions are summarized in Table 1. 

Table 1. 

Flame and Flow Conditions 


Central Jet Coflow 


Qh 2 (slm) 

25 

Qh 2 (slm) 

225 

Qn 2 ( slm) 

75 

Qair (slm) 

2100 

MajET 

0.3 


0.25 

T(K) 

300 

T (K) 

1035 

Vjet (m/s) 

101 

V coflow (m/s) 

3.5 

RejET 

23,600 

RecOFLOW 

18,600 

djET (mm) 

4.57 

Dcoflow (cm) 

21 


Simultaneous Raman-Rayleigh-LIF Measurements 

Multiscalar experiments were conducted in the Turbulent Diffusion Flame (TDF) 
laboratory at Sandia’s Combustion Research Facility. The temporally and spatially resolved 
simultaneous measurements of major species, minor species and temperature were made with 
Raman- Rayleigh scattering and laser induced fluorescence (LIF) techniques. The LIF systems 
measured the OF! and NO, while the Raman-Rayleigh system measured the N 2 , O 2 , H 2 O, H 2 and 
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temperature. Details on the experimental setup and the calibration techniques are presented 
elsewhere [Smith et al. 1995, Nguyen et al. 1995, 1996, Barlow et al. 1988, 1996, Nooren 1998], 

The system is setup to acquire data to produce joint statistics of the temperature and the 
major and minor species. The 5 lasers, (2 Raman-Rayleigh, 3 LIF) fire virtually simultaneously, 
with a delay of 100ns between shots. The laser pulse rate is 10 Hz. The spatial resolution of the 
system is 750pm. 

The precision and accuracy of the Raman-Rayleigh-LIF system is determined with the 
use of flat calibration flames [Barlow et al. 2000], The precision of single-shot measurements in 
a H 2 flame (no fluorescence interferences) is limited by the photoelectron shot noise [Dibble et 
al. 1987], The accuracy of the equipment is ultimately determined by the calibration flame 
measurements. Figure 2 shows the measurements of a CH 4 /Air Hencken-bumer taken prior to 
the H 2 /N 2 lifted flame experiment. The rms of the concentration or temperature and the mixture 
fraction data binds the ellipse that surrounds each point in Figure 2. 



Figure 2. 

Processed mean values of temperature and concentrations in the CH 4 -Air Hencken-burner 
flames. 


Probability Density Function Modelling Method 

Joint probability density function (PDF) methods have been demonstrated to be useful in 
the modeling of turbulent combustion with local extinction. Utilized is the joint probability 
density function for composition only. Equation (1) shows the instantaneous transport and 
production of N reactive scalars c^. 

+ P v * ^ = v * (p°V<k) + ,(^3 ( 1 ) 


The corresponding density-weighted PDF equation is: 
<P>§ + (P)2*VP + <p)| ^-(w,( Vi ,...., ¥n )p) 


= -V.| 


<P> 


9='lf ) p 


3 




( 2 ) 
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where \|/j are the sample space variables corresponding to the reactive scalars (J),; v and v" denote 
the density weighted average velocity and its variance. In Equation (3) £,, denotes the scalar- 
dissipation rate defined as. 

£ij = DVc^Vcl), (3) 

The turbulent flux and scalar dissipative terms appearing in the PDF transport equation are 
modeled. A gradient diffusion model is used for the turbulent flux and the Curl mixing model 
[Pope 1990] is used for the scalar dissipation rate. 

The Monte Carlo simulation technique is used to compute the transport equation for the 
pdf [Chen & Kollmann 1988], The Monte Carlo technique handles a collection of stochastic 
particles, involved in a simulation of convection, turbulent diffusion, molecular diffusion, and 
chemical reactions. The distribution of these stochastic particles in composition space then 
determines the form of the pdf. Subsequently, the single point statistics of the scalar variables (j), 
thru (|)n are determined via \)/i and the pdf. The computation power required by the multi- 
dimensional joint scalar pdfs limits the complexity of the reduced mechanisms for the 
combustion process. Currently, a 5-step reduced chemistry (4 steps for combustion and 1 step for 
thermal NO formation) is used for the modeling. In the reduced chemistry, both HO 2 and H 2 O 2 
are assumed to be in steady state, which may be relaxed in future modeling work. 

RESULTS AND DISCUSSION 


The structure of the H 2 /N 2 lifted jet flame is measured by obtaining several radial profiles 
and an axial profile of temperature and species concentration. Figure 3 shows the points where 
measurements were taken. The region in the flame where flame stabilization occurs was of 
interest; therefore several points were measured in this region. The axial profile consists of 
measurements from z/d=l to z/d=34 from the nozzle exit. Radial profiles were also obtained at 
several vertical distances (z/d = 1, 8, 9, 10, 11, 14 & 26). The radial domain covered by these 
radial sweeps was - 1 0mm to 50mm with step sizes typically between 1 and 3mm. 

The single shot data has been processed and the Favre averages and conditional mean 
averages along with their respective rms fluctuations were calculated. A modified Bilger 
formulation [Bilger et al. 1990], with the carbon components omitted since there is no 
appreciable carbon in the system, determines the mixture fraction. 


f = 

1 Bilger 


2M h ( Yh 
2M^( Yh >' “ 


Y »j-M;( Y o- Y °J 

Yh - 2 ) _ m;( y °> ~ Y °’ 2 ) 


( 4 ) 


Where Y is the elemental mass fraction and M is the atomic weight. Subscript 1 denotes the fuel 
stream and subscript 2 denotes the oxidizer stream. The elemental mass fractions at the 
boundary conditions are tabulated below in Table 2. The stoichiometric mixture fraction for 
these boundary conditions is fs = 0.475. 
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Stream 

Y h 

Yo 

Y n 

Fuel, 1 

0.0237 

0.0 

0.9763 

Coflow, 0 

0.0074 

0.2313 

0.7613 


Radial and Axial Profiles 

The radial profiles at z/d =1,14 and 26 of Favre averaged temperature and species mole 
fractions are plotted in Figure 4. The radial profiles at z/d =8, 9, 10 and 1 1 are included as an 
appendix. Two important design benefits are demonstrated in Figure 4. First, relatively flat 
profiles at z/d=l show that the boundary conditions are well defined. Second, the integrity of the 
coflow is undiminished as temperature and species measurements in the far field remain steady 
throughout the range of measurements in the vertical direction. 

Axial profiles of the Favre averaged temperature, mixture fraction and specie mole 
fractions are plotted in Figure 5. The results from the PDF model are also plotted with the 
experimental data. Two PDF model simulations were conducted, for both a reacting and non- 
reacting (pure mixing) flow. The solid lines indicate the reacting case and the dashed lines 
denote the non-reacting case. 


NASA/CR— 2002-2 1208 1 


6 








-B- XOH 
-O- XNO 




0 10 20 30 

r (mm) 




-B- XOH 
-O- XNO 




0 10 20 30 


r(mm) 


Figure 4. 

Radial profiles of Favre averages for temperature and specie concentrations at axial locations 
z/d=l (left column), z/d=14 (central column) and z/d=26 (right column). 


Lifted flames exhibit characteristics similar to both reactive and non-reactive flows, with 
the transition occurring at the base of the flame. The flame front acts as a barrier, limiting the 
mixing of the two streams. Therefore, as in the case of the mixture fraction, the non-reacting 
case exhibits steep gradients in the z-direction, in contrast to the more gradual trend in the 
reacting case. After the potential core is dissipated, the experimental results show this pure 
mixing and preheat in the region between the nozzle exit and the base of the flame. At the onset 
of combustion at the axis (z/d=13 ) the experimental results then demonstrate the limited mixing 
caused by the flame, the production of combustion products FLO, OH and NO and increase in 
temperature. However, the model predicts an accelerated ignition time, as can be seen in the 
temperature and OH mole fraction plots, this reduces the amount of mixing and entrainment of 
the coflow and central jet. 
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Axial profiles of the Favre average and rms fluctuation compared with the PDF model results. 
The connected circles are the experimental results. The reacting flow is plotted with solid lines 
while the non-reacting flow is plotted with dashed lines. The vertical dotted lines mark the onset 
of combustion at the axis ( z/d=13 ). 
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The hot oxidizer is entrained and preheats the fuel stream. The oxygen mass fraction 
reaches a maximum at approximately z/d=3 where the temperature is =550 K and the mixture 
fraction is =0.8. 

The aggregate experimental results of the radial profiles were used to generate a low- 
resolution color map of the temperature and species concentrations. The axial profile data was 
excluded from this process since it produced a significant bias for gradients in the radial 
direction. There is good resolution in the area where the flame stabilization occurs since several 
points were measured in this region. The temperature and OH color maps for the aggregate 
experimental data, and the results from the PDF model are plotted in Figure 6. Note the different 
YOH scales for the PDF and experimental results. 

25 
20 




Experimental and PDF color maps of temperature (top) and YOH (bottom). Note the different 


YOH scales for the PDF and experimental results. 
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Comparisons between the experimental and numerical results show that the numerical 
model does not accurately predict flame lift off. Sudden spatial increases in OH levels mark the 
flame edge [Tacke 1998], For our purposes here, we define the flame as the zone where the OH 
mole fraction is greater than 6 x 10" 4 . The experimental flame stabilization height is 
approximately H EX p/d=10 and the model predicts H PDF /d=2. Likewise, the radius of the base of 
the flame is Rs/d=1.5 and the model predicts RpDF/d=0.75. The levels of YOH in the experiment 
are approximately 50% of those predicted by the model. Above the stabilization point, the 
experiment and model agree. 

A reason for the difference in results between the experiment and the model is that the 
boundary conditions were slightly different. The raw experimental data goes through an iterative 
reduction process and preliminary experimental results were inadvertently used as boundary 
conditions for the PDF model. The differences between the experiment and model boundary 
conditions are listed in Table 3. These differences in the boundary conditions (2-3%) can have a 
substantial effect on the flame, therefore the PDF model must be modified and rerun. 


Table 3. 

Differences in preliminary and final experimental results 



Model 

Experiment 

Calculated* 

Jet 

Temperature (K) 

295 

291 

300 

xh 2 

0.254 

0.254 

0.252 

xn 2 

0.737 

0.746 

0.748 

Coflow 

Temperature (K) 

1070 

1035 

1080 

xo 2 

0.15 

0.147 

0.148 

xh 2 o 

0.102 

0.099 

0.102 

xn 2 

0.746 

0.753 

0.749 


* The calculated values are based on mass flow controller readings and equilibrium calculations. 

Conditional Statistics and Issues of Scalar Transport 

The central fuel jet entrains hot products from the coflow, thus evolving into a preheated, 
partially premixed flow prior to the onset of combustion. The Favre averaged measurements 
presented in the previous section show the entrainment of oxidizer into the fuel stream prior to 
combustion, thus further diverging the combustion statistics of this flame from those of a system 
at equilibrium. The scatter data of the temperature and oxygen mass fraction are presented in 
Figure 7. The data plotted in Figure 7 is a composite of the single shot data from the radial 
profile at z/d=14. Approximately 6000 data points are plotted. The partially premixed nature of 
the flame is shown by the penetration of oxygen to the rich mixture fraction space. 
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Figure 7. 

Scatter plot of single-shot temperature and YO 2 measurements at z/d = 14. The number of data 
points is approximately 6,000. The vertical dotted line indicates the stoichiometric value of the 
mixture fraction. 

The conditional means with the conditional rms (± 0 ) fluctuations plotted as error bars for 
the temperature and oxygen mass fraction corresponding to the scatter data in Figure 7 is plotted 
in Figure 8. The conditional means and rms are calculated at constant intervals of Af=0.02. 
Equilibrium calculation results are also plotted as solid lines. 



Figure 8. 

Conditional means and rms fluctuations plotted as uncertainty bars. The solid line is the results 
of equilibrium calculations of the fuel and oxidizer. 


Tsuji opposed flow laminar flame calculations were made with full molecular transport 
and equal molecular and thermal diffusivities. The results from the calculations with full 
molecular transport are plotted in Figure 9 along with the conditional means and rms fluctuations 
for z/d=14. 
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Fbilger 

Figure 9. 

Results of steady strained opposed-flow laminar flame calculations with full molecular transport. 
Conditional means and rms fluctuations ( error bars) of temperature, H 2 O and O 2 mass fractions 


are plotted against laminar flame calculations with strain rates of a=8,000 1/s (solid lines) and 


a=20,000 1/s (dashed lines). The vertical dotted line indicates the stoichiometric mixture 


fraction. 


There are considerable discrepancies between the full molecular transport calculations 
and the experimental results. This is in agreement with the well-documented findings of Barlow 
et all 996 & 2000. The mixture fraction domain over which reactions occur is smaller than what 
the calculations predict. This shows that the turbulent mode of mixing may be dominant. 
Therefore, the calculations were rerun with equal diffusivities. 

The presence of oxidizer in the fuel stream is evident in the mixing line that diverges 
from the equilibrium curve in the rich mixture fraction space (Figure 8). Since the flame is 
lifted, the fuel stream becomes partially premixed, thus changing the effective boundary 
conditions for the opposed flow flames calculations. For the equal diffusivity calculations, the 
fuel side boundary condition was set to f ~ 0.8 which is approximately the mixture fraction 
corresponding to the onset of combustion as presented in the axial profiles (z/d=13). The 
temperature is approximated to be 440K, the temperature corresponding to f=0.8 on the pure- 
mixing line in Figure 8. The results of these calculations are plotted in Figure 10. 
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Figure 10. 

Results of steady strained opposed-flow laminar flame calculations with equal specie and 
thermal diffusivities. Conditional means and rms fluctuations (error bars) of temperature, H 2 0 
and 0 2 mass fractions are plotted against laminar flame calculations with strain rates of a=l,000 
1/s (solid lines) and a=4,000 1/s (dashed lines). The vertical dotted line indicates the 
stoichiometric mixture fraction. 

There is good agreement between the experimental and equal diffusion calculations. The 
contrast of this agreement with that of the experimental and full transport calculations suggests 
that turbulent transport is dominant over molecular transport. There is some discrepancy on the 
rich mixture fraction space. This may be due to the heat loss to the combustor in the opposed 
flow flame calculations. Future calculations will be performed to see if better agreement can be 
achieved by a more adiabatic system with increased flow velocity with constant strain rate. 

CONCLUSIONS: 


Simultaneous multiscalar measurements of a lifted H 2 /N 2 jet have been presented and 
provide the basis for insight into flame structure of a lifted flame. The flame exhibits the partial 
premixing inherent in lifted-off flames, while also being preheated by the entrained hot coflow. 
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Laminar flow calculations with full molecular transport have also been made and do not agree 
well with the experimental results. However, the laminar flame calculations with equal 
diffusivities agree well when the preheating and premixing in the jet prior to the flame 
stabilization is accounted for in the boundary conditions. An exploratory PDF simulation of the 
system has been conducted and has been presented. It does model the axial profiles somewhat 
accurately, but it does not model the lift-off height. 
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